Finite-size behavior near the critical point of QCD phase-transition 
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It is pointed out that finite-size effect is not negligible in locating critical point of QCD phase 
transition at current relativistic heavy ion collisions. The finite-size scaling form of critical related 
observable is suggested. Its fixed point behavior at critical incident energy can be served as a 
reliable identification of critical point and nearby boundary of QCD phase transition. How to find 
experimentally the fixed point behavior is demonstrated by using 3D-Ising model as an example. 
The validity of the method at finite detector acceptances at RHIC is discussed. 
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Quantum colordynamics (QCD) has predicted quark 
deconfinement and chiral symmetry restoration at finite 
temperature and density Lattice-QCD has shown 
that the transition is crossover at vanishing baryon chem- 
ical potential /xb 0- The QCD based model indicates 
that the crossover turns to be a first-order phase tran- 
sition at larger values of fiB H- The endpoint of the 
first order phase transition to the crossover is referred as 
critical endpoint, or critical point. All these show that 
the transition from hadron phase to quark-gluon-plasma 
(QGP) phase can happen in 3 possible ways. 

The data from current relativistic heavy ion experi- 
ments show that the QGP has been formed at RHIC en- 
ergies [|[. But the position of critical point in QCD phase 
diagram is not clear from the theoretical side. Well de- 
fined character of critical point, the divergence of corre- 
lation length, warrants the possibility of finding it exper- 
imentally. The goal of the beam energy scan at RHIC, 
and future heavy ion experiment at FAIR is aimed to 
pass through the critical incident energy. 

In relativistic heavy collision, both the size and dura- 
tion of formed system are finite. From event to event, the 
overlapping area, i.e., the formed system size, varies with 
impact parameter. For finite-size system, critical behav- 
ior changes with system size (L). If the system size is too 
small, the correlation length can not be fully developed 
to cause a phase transition. If the system size is large 
enough and the correlation length (£) is much smaller 
in comparison with system size, the system can still be 
considered as infinite large. The critical behavior under 
thermal limit is available. This is why the non-monotonic 
behavior is suggested as an indicator of critical point 0- 
0] in a long period. 

Nevertheless, non-monotonic behavior is not unique to 
critical point. In case of first order phase transition, or 
crossover, some observables also show the non-monotonic 
behavior Q. The absence of non-monotonic behavior 
does not exclude the existence of critical point, such as 
the maximum cluster size in 3D-Ising model shown in 
Fig. 1(a). 

Moreover, if the correlation length is comparable with 
system size, the finite-size effect is not negligible. When 



correlation length is large than | of system size, it has 
been shown that the finite-size effect has to be taken into 
account [1, 0. 

Although, it is still difficult to estimate the size of the 
formed system and correlation length at critical point 
in relativistic heavy ion collisions. A rough estimation 
shows that the system size at freeze-out is less than 
12fm [13, [TTj] - The correlation length is round 6fm for 
typical nuclear collsions [TJ E3- After considering the 
finite evolution time, or finite-size, it is argued that the 
maximum of correlation length may not be beyond 2-3fm 
at critical point [13hT5| . Base on those estimations, the 
ratio of correlation length to system size is round 
which is in the region larger than i. So in current rela- 
tivistic heavy ion collisions at RHIC, the finite-size effect 
most probably has to be taken into account, rather than 
negligible 
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In accounting for the finite-size effects, the critical be- 
havior of all suggested observables, such as the fluctu- 
ations and correlations of transverse momentum, mul- 
tiplicity, conserved charges [lH, and in particular, the 
higher order moments of baryon numbers [171 ] , should be 
re-examined under the frame of finite-size behavior of 
critical point and nearby boundary. 

In the letter, we first discuss the finite-size behavior 
of the critical point, the first order phase transition, and 
the crossover in general. Secondly, we suggest the finite- 
size scaling form of the critical related observable in rel- 
ativistic heavy ion collisions. Its fixed-point behavior at 
critical incident energy can be served as a reliable iden- 
tification of critical point and nearby boundary of QCD 
phase transition. Thirdly, we demonstrate how to locate 
the fixed point from experimental observable by using 
3D-Ising model as an example. Finally, the validity of 
the method at finite detector acceptance at RHIC is dis- 
cussed. 

For the second order phase transition, the critical 
behavior is well described by finite-size scaling. It 



was firstly proposed from phenomenological [181 ] and 
renormalization-group [l9[ theories, and was approved by 
the Monte Carlo results of finite systems in different uni- 
versal classes [2(J. This scaling form not only describes 
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the behavior of the observables at different system sizes, 
but also indicates the position of the critical point and 
the critical exponents in infinite system. Therefore, from 
the finite-size scaling of critical related observables, the 
position and critical exponents of critical point can be 
precisely extracted. This has been implemented in locat- 
ing critical point of multi-fragmentation nuclear liquid- 
gas phase transition (2l| . 

In contrast to the critical point, the finite-size behavior 
of first order phase transition has not been well under- 
stood in general [22j ■ But the finite-size scaling behavior 
of first order phase transition is shown to correspond to 
so-called discontinuity fixed points of the renormalization 
group transformations, which are characterized byeigen- 
value exponents equal to the spacial dimension [23j. So 
the finite-size scaling form pertains, and the scaling expo- 
nents are the spatial dimension, in contrary to the critical 
exponents of critical point. The phenomenological theory 
of finite-size scaling at first-order phase transition is pro- 
posed by K. Binder and D.P. Landau and found to be in 
good agreement with Monte Carlo simulation results [24| . 

Different from the critical point and the first order 
phase transition, at the crossover region, there is no sin- 
gularity in all kinds of observables. The observables are 
system size independent @, [25|. But it should be no- 
ticed that this holds only when the system size is not 
too small. When the system size is very small and the 
finite correlation length is comparable with the system 
size, the observables will become larger and larger when 
the system size goes to smaller and smaller. 

In heavy ion collisions, the critical related observables 
are generally considered to be the fluctuations of con- 
served charges, like baryon number, electric charge, and 
strangeness [III E3, IHj]. The incident energy and cen- 
trality dependence of some related observables are fully 
investigated in current heavy ion experiments (27j . 

Incident energy y/s is a controlling parameter, like tem- 
perature T, or external field h in thermodynamic sys- 
tems. The centrality, i.e., impact parameter, presents 
the overlapped area of two incident nuclei. It directly 
related to the size of the formed system, and randomly 
fluctuates from event to event. 

So the finite-size scaling in nuclear collisions can be 
generalized as following. When the size of the formed 
matter L is much larger than the microscopic length 
scale (which is less than f fm) and incident energy is near 
the critical one v / s c , the critical related observable, e.g., 
Q[\fs , L) in general, can be written in a finite-size scaling 
form [18l - f20j ] . 

Q(V~s,L) = L x ^F Q (rL 1 /n- (1) 

Where r = (\/s — \/^ c )/\/^c * s reduced incident energy, v 
and A are the critical exponents of the correlation length 
£ = £o T v an d the observable, respectively. They charac- 
terize the universal class of the phase transition. Finite- 
size scaling indicates that the observable at different sys- 
tem sizes can be re-scaled to an identical scaling function 
Fq with scaled variable tL x I v . 
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FIG. 1: (Color online) Upper panel: (a) and (b) are the tem- 
perature dependence of maximum cluster size and maximum 
cluster size scaled by L d ~ 13 ^ 1 ' , and (c) is the scaling function 
of maximum cluster size in 3-D Ising model at different lattice 
sizes L. Lower panel: the same measures as in corresponding 
upper panel, but in a sub-system containing fifty percent of 
the lattice sites. 



At critical energy, y^s = v^ C ' tne sca l e d variable 
{tL 1 / v — 0) is independent of system size L, and the 
scaling function becomes a constant, 

Fq(0)=Q(VS c , £)£"*/". (2) 

It shows that the fluctuation of critical related observ- 
able is self-similar at different size scales. In the case, 
the energy dependence of the observable at various of 
system sizes will intersect to this point, i.e., fixed point. 
The energy of the fixed point indicates the critical inci- 
dent energy. As an example, we show in Fig. 1(b) the 
fixed point behavior of the maximum cluster in 3-D Ising 
model, which is su ppo sed to be the same university of 
the deconfincment [28|. We can see that the maximum 
cluster at different lattice sizes intersect exactly at the 
fixed point, i.e., critical point. 

Reversely, if we can find the fixed point from incident 
energy dependence of properly scaled observable in heavy 
ion collisions, which are measured at different centrali- 
ties, i.e., system sizes, it will indicate the existence of 
critical point. 

In order to find the exponent of the scale, and incident 
energy of the fixed point, we can first present the incident 
energy dependence of critical related observable at differ- 
ent system sizes, similar to Fig. 1(a). Then multiply a 
size factor to the observable Q{-\fs, L), i.e., Q(y/s, L)L~ a , 
and change the parameter a from — oo to oo to see if all 
size curves interact to a point for a certain value of ao at 
a certain incident energy, e.g., Fig. 1(b). 

In experiment, the point liked behavior can be quan- 
tified by the width of all size points. At a given incident 
energy, the width is usually defined as the square root of 
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FIG. 2: (Color online) Incident energy dependence of point 
liked behavior, -D m in(\/s 5 a o) vs. y^s. 



X 2 of all size points, i.e., 



N L -1 



(3) 



Nl is the number of points, and Xq/^ jj^-a i s error 
weighted variance of all size points, 
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[Q(^,Li)LT a -(Q(^,L)L- 
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w,i = 5 [Q(\/~s, Li)L^ a ~\ are the experimental error of 
[Q(y r s,Li)LY a ], where both the errors of the observ- 
able Q(yfs,Li) and system size LJ a contribute to. 
(Q(y/s, L)L~ a ) is also error weighted mean, 



{Q{yfs,L)L- a ) 



Ef=r iK 
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For example, in fig. 1(b), this width at a given tempera- 
ture is the distance between two violet arrows. Therefore, 
if at a given incident energy, the minimum of D{^/s : a) 
is round 1 at a , i.e., D m i D (^/s, a ) ~ 1, it can be recog- 
nized as an experimental point. While, if it keeps larger 
than 1, there is no point liked behavior. 

For QGP formed system [J], the following 3 cases 
should be expected. (1) D m i D (y/s, ao) at a certain inci- 
dent energy is round 1, and at nearby incident energies, 
it is always larger than 1, and corresponding ao is not 
an integer, as green curve shown in the middle of Fig. 2. 
This may indicate the existence of the fixed point, i.e., 
the critical point in Eq. ^j. The critical incident energy 
is the energy of the fixed point and the obtained param- 
eter ao is the ratio of critical exponents, i.e., \/v = ao, 
cf. Fig. 1(b). 

In the case, the critical behavior should be further con- 
firmed by the scaling function, 



F Q (rL 1 /»)=L- a °Q(^,L). 
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Here the critical exponent of correlation length v is a 
fitting parameter. If the data at all incident energies and 



system sizes can be well fitted by the scaling function, 
the critical point and the critical exponents are finally 
determined, cf., fig. 1(c). 

(2) D m i n (y/s, ao) at a certain incident energy is round 
1, and at nearby incident energies, it is always larger 
than 1, and corresponding ao is an integer as green curve 
shown in the left side of Fig. 2. This indicates also the 
existence of the fixed point, but scaled power is trivial 
integer. It implies the region of the first order phase 
transition. The incident energy of the fixed point is the 
transition energy of the first order phase transition. The 
scaling function of the observable should be simply for- 
mulated by the spatial dimension, instead of the critical 
exponents in Eq. (1). 

If a is zero, there are two possibilities. It could be the 
critical point with critical exponent A = 0, like Binder 
cumulant ratio 29] , or the region of the first order phase 
transition. The final identification is their specified scal- 
ing functions, as discussed above. 

(3) D m i n (y/s, ao) is round 1 at all incident energies, as 
green line shown in the right side of Fig. 2, and corre- 
sponding ao is an integer. This indicates all size curves 
are overlapped, and there is in fact no fixed point. It 
corresponds to the transition of crossover. 

It should be stressed that the observables we men- 
tioned here are the intensive variable, like susceptibility. 
If the observables are extensive variable, such as the fluc- 
tuation of the particle number, ((N — N) 2 ) = TVx, the 
trivial size dependence are included, and can be merged 
to the power a. 

The size of the formed matter in heavy collisions is 
mainly determined by overlapping area of two incident 
nuclei. This area is proportional to the number of par- 
ticipant nucleons and is quantified as ccntrality. The ini- 
tial size of the formed matter can be approximately esti- 
mated by the square root of the number of participants, 
\J -^Vpart ■ The maximum size is ^/2Na, Na is the number 
of nucleons of incident nucleus. The ratio, 
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presents the relative size of the initial system. 

The system size V at transition should be larger than 
initial one L and monotonically increase function of L, 
i.e., L' = cL 1+s with S > in general. Whether we take 
V or L in Eq. ([T]), the scaling exponents will be different, 
but the position of critical point will be the same. So 
the initial size is a good approximation in locating the 
position of critical point. 

It should also be noticed that the detectors at cur- 
rent relativistic heavy ion experiments cover a part of 
the phase space, and only a part of final state particles 
is accepted. Even if the critical related information are 
survived in the final state observables, whether the finite- 
size behavior of detected subsystem is preserved has to 
be examined further. 

The finite size behavior of a sub-systems is demon- 
strated in 3-D Ising model. The size of sub-system is 
chosen to be a certain percent of the whole lattice sites. 
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Changing the lattice of the whole system, the effective 
sites of the sub-system, L e ff, vary with it. We find that 
the finite size behavior of sub-system keeps valid as long 
as the size of sub-system is within the range of finite size 
scaling. 

In the lower panel of Fig. 1, the finite size behavior of 
the maximum cluster size at various £ e // are presented. 
Where the size of sub-system is 50% of the whole sys- 
tem. In comparison with the corresponding results of the 
whole system shown in upper panel of Fig. 1, the suscep- 
tibilities of the sub-system is different from that of the 
whole one, but the position of fixed point indicates the 
same critical temperature, T c = 4.51J. Moreover, the 
maximum cluster at different sub-system sizes are well 
scaled to an identical scaling function. Therefore, the 
suggested finite size behavior should be visible at a de- 
tector with a relative large acceptance, like RHIC/STAR. 

In the summary, we point out that the finite-size ef- 
fects are not negligible in locating the critical point of 
QCD phase transition at current relativistic heavy ion 
collisions. At the crossover, critical point and first order 
QCD phase transition, the finite-size scaling behaviors of 
the critical related observable are suggested. 

The critical point of QCD phase transition can be 
found by the appearance of the fixed point with a non- 
integer power in scaled size factor, and the finite-size scal- 



ing function of the observable. The region of the first or- 
der phase transition is identified by the fixed point with 
an integer power in scaled size factor and the scaling 
function which is determined by spatial dimension. 

At the region of the crossover, the behavior of the fixed 
point is absent, and the scaling function reduces to the 
incident energy dependence of observable, which is sys- 
tem size independent. 

At a given incident energy, the width of the observ- 
ables at various centralities is suggested as a quantifi- 
cation of point liked behavior. The energy dependence 
of the width at different orders of phase transitions are 
shown. When incident energy scans from high to low, the 
deviation of minimum width from point like behavior will 
indicate the appearance of the critical point. 

Finally, for a finite acceptance detector, we demon- 
strate that the finite-size behavior of critical related ob- 
servables keep valid as long as the detected subsystem is 
large enough. 
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